Journal of Biomechanics

Background and ObjectiveThe foot-ankle complex is a highly articulated and mechanically constrained system, often simplified as a chain of few rigid segments, neglecting many bone-to-bone motions and raising questions about the accurate representation of interaction with ground. This study proposes a new reduced-order multibody formulation that captures intrinsic kinematic constraints of the foot through motion synergies. MethodsBones kinematic coupling, or motion synergies, were experimentally derived from weight-bearing CT scans using principal component analysis. These couplings were embedded in a synergy-based multibody kinematic optimization framework describing the foot-ankle with five degrees of freedom: ankle flexion; foot adduction, pronation, and arching; and toe flexion. Model accuracy was evaluated against bone-level experimental kinematics. The model was applied to gait data from healthy, flat, and diabetic feet and compared with a standard multi-segment foot model, assessing robustness by progressively reducing the number of skin markers. ResultsAverage errors were about 1{degrees} and 0.5 mm when using subject-specific synergies and below 7{degrees} and 4 mm when scaling the generic model, matching or exceeding the accuracy of existing models. Reliable reconstruction was obtained using only four foot markers. In clinical gait analysis, the model showed superior discrimination between populations and enabled assessment of transverse arch deformation, not accessible with conventional models. ConclusionThe proposed synergy-based model provides an accurate, low-complexity framework for reconstructing bone-level foot and ankle kinematics, substantially simplifying gait analysis while improving biomechanical interpretability. This framework supports future integration with dynamic models aimed at studying load transmission in the foot.

Importance: Foot drop impairs mobility for many children globally, causing life-long health issues. Existing treatments are costly, custom-made, and require frequent clinical visits. A new, low-cost, off-the-shelf splint (OrthoPed) could improve access and user experience. Objective: To determine the feasibility of recruiting children (4-17 years) with moderate foot drop and collecting biomechanical, clinical, and patient-reported outcomes to compare OrthoPed with existing treatments. Design: Single-centre cross-sectional feasibility and pilot study informing a future randomised clinical trial. Participants: Twelve children (target=20; mean age=10.6 {+/-} 3.5 years; 2 females) with moderate foot drop and prescribed orthotic support were recruited via physiotherapy. Intervention: The new OrthoPed splint was compared against existing treatments: ankle foot orthoses (AFOs) and Lycra socks. Main outcome measures: Primary outcome: recruitment and retention rates. Secondary outcomes: biomechanical and clinical gait measures, alongside useability and performance questionnaires. Results: Recruitment reached 22% of eligible participants (an "amber" rating for future trials). Despite four dropouts due to treatment burden, all outcome measures were successfully collected. Preliminarily, OrthoPed supported more natural gait mechanics than AFOs and offered better usability and comfort than AFOs and Lycra socks, potentially enhancing adherence. Conclusions: Recruiting children for orthotic trials is feasible, though coordinating gait testing with routine clinical appointments could improve future recruitment. Importantly, low-cost orthotic devices may provide better usability, accessibility and adherence than existing prescribed options.

Conventional diagonal stride skiing traditionally includes a glide phase, characterised by a period of relatively passive gliding on one ski. While the glide phase may take advantage of low ski-snow friction, it does not exhibit the same whole-cycle mechanical energy fluctuations seen in running or walking on foot. A new sub-technique, known as running style, substantially reduces the glide phase and may alter the role of elastic tissues, making the movement pattern more similar to uphill running on foot in its temporal organisation. We examined knee extensor and plantar flexor muscle-tendon behaviour in eight competitive skiers performing conventional diagonal and running techniques on a treadmill inclined at 10{degrees}. Using synchronised ultrasonography, 3D kinematics, ski forces and EMG, we quantified gastrocnemius medialis and vastus lateralis fascicle and muscle-tendon unit (MTU) dynamics in both the running (RUN) and conventional (CON) styles. Shorter glide and total cycle durations during RUN shifted MTU peak length and velocity earlier during the kick phase. Fascicles in both muscles operated at similar velocities across techniques, showing MTU-fascicle decoupling. Vastus lateralis fascicles shortened at higher absolute peak velocities than gastrocnemius in both conditions, while normalised velocities were similar. RUN increased preactivation and advanced EMG timing, while integrated EMG during the kick was lower compared to CON. These findings suggest that, despite large shifts in external mechanics between glide-based and more running-like skiing, elastic tissues may help stabilise fascicle behaviour and preserve a similar contractile strategy across muscles and techniques.

Background: Quadriceps strength and landing mechanics are two modifiable factors associated with anterior cruciate ligament (ACL) injury risk. Collecting detailed biomechanical data is an arduous task. Identifying a relationship using more easily measured variables, such as quadriceps strength, would offer value for athlete counseling and injury prevention programs. Although quadriceps weakness has been associated with altered landing strategies in ACL-reconstructed (ACLR) individuals, this relationship is less clear in healthy athletes. Purpose: To investigate the association between isokinetic quadriceps strength and peak knee flexion angle during a vertical drop jump in healthy adolescent athletes. Study Design: Secondary analysis of previously collected data. Methods: Healthy adolescent athletes had their dominant leg quadriceps strength measured using an isokinetic dynamometer at 60{degrees}/s from 0-90{degrees} of knee flexion. Landing mechanics were assessed during a vertical drop jump using three-dimensional motion capture synchronized with force plates. Pearson correlation was used to evaluate the association between quadriceps strength and peak knee flexion angle during landing, with statistical significance defined as p < .05. Results: There was a weak negative correlation between quadriceps strength and peak knee flexion angle (p = .017, R = -.22 [-.04, -.38]), suggesting that stronger athletes achieved greater knee flexion angles. Discussion: Greater quadriceps strength was associated with increased peak knee flexion angles during landing; however, the weak correlation suggests that strength explains only a small portion of the variability in landing mechanics. These findings deviate slightly from prior literature in healthy populations but are consistent with studies demonstrating that greater quadriceps strength is associated with achieving greater peak knee flexion in ACLR patients. Accordingly, quadriceps strengthening should remain a key component of multifactorial ACL injury prevention programs.

Gait asymmetry is a common manifestation of walking impairment among clinical populations. We recently developed a novel treadmill walking approach called dynamic treadmill walking that can provide asymmetric gait training by changing the treadmill speed between fast and slow speeds within a single stride. Here, we studied the energy expenditure associated with a variety of dynamic treadmill walking conditions. We hypothesized that the metabolic power required for dynamic treadmill walking in all conditions would approximate the metabolic power associated with conventional walking at the mean of the fast and slow speeds employed in the task. Eleven young adults without gait impairment walked on an instrumented treadmill and breathed into a metabolic measurement system. During dynamic treadmill walking, the treadmill fluctuated between 0.75m/s and 1.50m/s, each for 50% of an individuals stride time. We used a metronome to synchronize participants right heel-strikes with four different timing conditions. Net metabolic power during dynamic treadmill walking was significantly greater than normal walking at the mean speed of the task (1.125m/s) and generally lower than walking at the fast speed (1.5m/s). We did not observe any significant associations between net metabolic power and several measures of gait asymmetry during dynamic treadmill walking. These findings establish dynamic treadmill walking as a promising technique for improving gait symmetry in individuals who cannot tolerate fast treadmill walking, a common gait rehabilitation approach. Future work will assess the feasibility, metabolic demands, and clinical efficacy of using dynamic treadmill walking to improve gait symmetry in clinical populations. Key Points SummaryO_LIDynamic treadmill walking (i.e., walking with oscillating treadmill speeds) has previously been shown to drive gait asymmetries, but little is known about the energy expenditure required to complete the task. C_LIO_LIOur hypothesis was that dynamic treadmill walking would have similar metabolic power requirements to normal walking at a speed that is intermediate between the two dynamic treadmill walking speeds. C_LIO_LIWe found that dynamic treadmill walking actually requires metabolic power that is greater than the average of the two belt speeds, but less than that used for fast walking. C_LIO_LIDynamic treadmill walking is a promising and clinically translatable technique for rehabilitating populations with gait asymmetries that is not more energetically costly than fast treadmill walking, a common gait rehabilitation approach. C_LI

Biomechanical simulations of complex hand motions remain scarce, due to challenges that span computation and data acquisition. Using a computer vision-based motion capture approach, a 23-degree of freedom musculoskeletal model, and direct collocation optimization, we performed muscle-driven simulations to track hand kinematics from 7 participants performing American Sign Language gestures. While proximal joints were tracked accurately, interphalangeal joint tracking was significantly worse, with a consistent flexion bias. Modifications to finger extensor muscle paths that incorporated the dual-inserting nature of the extensors improved accuracy, suggesting better representation of extensor force distribution across distal joints may be necessary for accurate hand simulations.

Aortic dissection is life-threatening due to continued loss of medial integrity that may culminate in secondary rupture within hours to days. While pre-existing defects or hemodynamic loads compound structural deterioration of the aorta, pathological progression from symptomatic dissection channel to lethal transmural tear is poorly understood. We examined the structure of referent and acutely dissected ascending aortas by microscopy. Elastic, collagen, and cellular components of non-dissected media were intricately interconnected. Medial damage in dissection lesions was traced from ingress to central to peripheral areas. Entry tears broke cleanly through successive laminae leading to cavernous false lumens in which medial structure was destroyed. Nearby laminae with widening between flanking elastic lamellae (termed minor delaminations) were filled with blood and showed severe medial damage. Farther laminae without delamination but containing red blood cells (termed blood extravasation) displayed moderate medial damage. More distant, non-delaminated laminae with accumulation of albumin but not red blood cells (termed plasma extravasation) exhibited mild medial damage. Varying medial hemorrhage with scattered sloughing of laminae was observed along the entire false lumen. We conclude that hydraulic fracturing of residual dissected media by pressurized blood via communications from the false lumen contributes to further structural weakening of the aortic wall.

BackgroundThe pivot shift (PS) test is the most specific clinical examination for anterolateral rotational instability in ACL-deficient knees, yet grading remains subjective, as evidenced by poor inter-observer reliability, particularly for Grade 2. Since low-grade (Grade 1) versus high-grade (Grades 2/3) PS is the threshold for recommending lateral extra-articular augmentation, performing the test in awake clinic patients limits grading reproducibility and introduces variability in surgical decision-making. Existing methods to quantify the pivot shift usually require examiner-performed testing under general anaesthesia. No prior approach has ascertained PS grading from a separate patient-performed functional movement. PurposeTo evaluate the feasibility of a machine learning (ML) classifier, trained on kinematic ultrasound bone-tracking signals acquired during a patients sit-stand-sit (SSS) knee movement, to predict their PS grade, and to clinically validate its ability to differentiate low versus high-grade PS. MethodsUltrasound bone-tracking kinematic data were collected during SSS manoeuvres in 23 ACL-injured patients using the GATOR device, and ground truth PS grades (0-3) were assigned under general anaesthesia by fellowship-trained orthopaedic sports surgeons. From the data collected, Leave-one-out cross-validation (LOOCV) was used to train the ML classifier. Clinical SSS data from 6 ACL-deficient patients was used for independent held-out validation of their low-grade (Grade 1) versus high-grade (Grade 2/3) PS. Multiple deep learning architectures (XceptionTime, InceptionTime, FCN, ResNet, ResCNN) and training strategies (including mixup augmentation and supervised contrastive learning) were tested. Performance was measured by one-versus-rest (OVR) AUC under LOOCV and by AUC (low vs high grade PS) from the held-out patient sessions. ResultsThe ML classifier achieved a maximum OVR AUC of 0.928 {+/-} 0.084 under LOOCV. Classifier performance increased with pivot-shift severity: Grade 3 was identified most reliably (AUC ~0.81; sensitivity 0.70-0.80), whereas Grade 2 remained the most challenging boundary (sensitivity 0.20-0.75 across configurations). For the clinically relevant binary classification of low-versus high-grade pivot shift, the classifier generalised well to a completely unseen patient cohort (AUC 0.889; accuracy 0.860; sensitivity 0.850; minimum-class sensitivity 0.767). ConclusionThe study demonstrates that kinematic ultrasound bone-tracking during sit-stand-sit contains transferable information about rotational instability severity in ACL-deficient patients, and represents the first reported approach to predict pivot shift grade from a patient-performed functional movement. The strong cross-validation performance confirms that the signals contain meaningful PS grade-discriminative information, but larger datasets targeting 50-100 sessions per grade will be required to achieve patient-level generalisation and advance this novel rotational instability assessment tool toward full clinical adoption. Level of EvidenceLevel IV, diagnostic feasibility study.

Various ankle-foot conditions (e.g., fractures, diabetic foot ulcers, and post-surgical recovery) require periods of complete non-weightbearing followed by gradually increasing partial loadings. However, existing assistive devices often provide inconsistent or uncomfortable offloading during gait. Additionally, prolonged proximal leg offloading can contribute to muscle atrophy, reduced bone density, and overuse of other body segments. We present a novel offloading ankle-foot orthosis (OLAFO) designed to overcome these limitations. The OLAFO features a patient-specific load-bearing shank brace, designed through a digital workflow and fabricated from a 3D-printed core reinforced with carbon-fiber composite lamination. Interlocking serrated side struts, adjustable in 2 mm increments, modulate load sharing between the shank and plantar surfaces. Furthermore, the OLAFO incorporates contact plates with a rocker profile informed by roll-over-shape measurements to support forward progression and gait symmetry. Proof-of-concept biomechanical verification in one able-bodied participant evaluated complete offloading, five partial-loading levels, and normal gait using a pressure walkway to compute vertical ground reaction forces and impulses. In complete offloading, the affected foot generated no contact pressures. Across partial-loading levels, the foot impulse increased from 14% to 53% of the total load and scaled linearly with strut height adjustments, supporting clinician-prescribed loading increments. Contralateral stance duration increased only modestly compared to commonly used assistive devices, indicating reduced compensatory loading on the intact limb. These findings demonstrate the proof-of-concept feasibility of the OLAFO, highlighting its potential for verifying full offloading and prescribing partial-loading targets during rehabilitation. Future research will evaluate performance across patient populations and clinical rehabilitation tasks.

ObjectivesTo determine whether early functional severity at presentation explains variability in return to sport (RTS) after ankle sprain in young athletes, compared with sprain subtype and injury mechanism. DesignRetrospective cohort study. MethodsAthletes aged [&le;]22 years with acute ankle sprains were identified from a prospectively maintained institutional database. Surgically treated cases were excluded. Functional severity at presentation was classified into three grades based on the ability to continue sports participation and ambulate immediately after injury. Injury mechanisms were categorized as high-energy deceleration (HED) or non-HED. RTS was analyzed as time to return and as prolonged RTS ([&ge;]4 weeks). Multivariable logistic regression was performed to identify factors independently associated with prolonged RTS. ResultsA total of 437 cases were included. Median RTS was 2.0 weeks (interquartile range, 0.0-4.0), and prolonged RTS occurred in 33.0% of cases. RTS duration increased stepwise with greater functional severity (p < 0.001). In multivariable analysis, functional severity was strongly associated with prolonged RTS (Grade 2: adjusted odds ratio [OR], 3.58; 95% confidence interval [CI], 2.07-6.19; Grade 3: adjusted OR, 24.53; 95% CI, 10.67-56.43; p < 0.001), and age was also independently associated (adjusted OR, 1.19 per year; 95% CI, 1.11-1.27; p < 0.001). Sprain subtype and injury mechanism were not independently associated with RTS after adjustment. ConclusionsEarly functional severity at presentation is the primary determinant of RTS after ankle sprain in young athletes. Apparent differences related to sprain subtype and injury mechanism are largely explained by initial functional impairment.

Thin structures such as heart valves and aortic dissection flaps interact dynamically with blood flow in human vessels. Their flexibility and capacity for large deformations generate complex, highly transient hemodynamic patterns over the cardiac cycle. Accurately resolving these interactions remains challenging for conventional boundary-fitted fluid-structure interaction approaches. We present an immersed boundary method for simulating thin structures in incompressible flow on unstructured grids. The method couples a stabilized finite element fluid solver with a nonlinear, rotation-free shell formulation through a direct forcing immersed boundary approach. The framework supports both weak (explicit) and strong (implicit) time-coupling strategies, enabling stable simulations over a wide range of solid-to-fluid density ratios. Hydrodynamic forces acting on thin structures are computed from fluid solutions sampled on both sides of the structure, allowing accurate force reconstruction for zero-thickness shells. To our knowledge, this is the first immersed boundary formulation that couples an unstructured finite element fluid solver with a two-dimensional, rotation-free shell model to simulate interactions between thin structures and incompressible flow. Fluid-structure coupling is achieved using predefined finite element shape functions, which provide consistent projection between Eulerian and Lagrangian fields without additional interpolation procedures. The framework is validated using three-dimensional benchmark problems involving thin structures. Then, valve-like model is used to compare strong and weak coupling strategies. Finally, the method is applied to an idealized type-B aortic dissection model. The proposed approach is implemented within the open-source software CRIMSON, a finite element platform for cardiovascular simulation.

ObjectiveRecently, alternatives to animal testing, such as new approach methodologies, are being developed in the orthopedic research field; animal models still provide valuable insights into the pathogenesis of knee osteoarthritis (OA). However, commonly used models develop OA much more rapidly and severely than those observed in human patients. We aimed to develop a novel murine model that closely mimics the slow progression of human OA with posterior Cruciate ligament (PCL) rupture. Design12-week-old C57BL/6 mice were induced to PCL-rupture (PCL-R) by manually applying an external tibial posterior translation force. We analyzed joint kinematics, histological observations, and bone structure to confirm the absence of concurrent injury on day 0. Then, joint stability and the pathophysiological progression of knee OA were analyzed at 8, 16, and 34 weeks post-PCL-R. The destabilized medial meniscus (DMM) model was also analyzed to compare the OA progression. ResultsNon-invasive PCL-R intervention induced the complete rupture in the central region of PCL without concurrent injury. The PCL-R group showed larger posterior tibial deviation than the INTACT (P=0.008). Regarding the range of motion in the PCL-R group, there was no limitation in range of motion on day 0, but extension limitations occurred at weeks 16 and 34 weeks. Histologically, articular cartilage degeneration in PCL-R was milder than DMM. In the subchondral bone, micro-CT reconstruction images indicated that, compared with the INTACT group, the DMM group observed progressive subchondral bone formation from 16 weeks post-surgery. In contrast, the PCLR group maintained the subchondral bone structure even at 34 weeks. ConclusionsPCL-R model induced mild abnormal mechanical stress depending on posterior instability, and cartilage degeneration occurred more slowly in this model than in DMM models.

Understanding how biological shape and movement interact with surrounding fluids represents a fundamental challenge at the intersection of biology, physics, and engineering. Fish locomotion exemplifies this challenge: body morphology and swimming kinematics together determine the hydrodynamic forces and flow structures that enable efficient propulsion and maneuverability. Whereas biologists have long sought to connect morphological variation to swimming performance, traditional morphometric approaches provide limited insight into the fluid mechanical consequences of shape differences. Similarly, although computational fluid dynamics can reveal detailed flow physics, simulating hydrodynamics across diverse and dynamic morphologies remains prohibitively expensive for systematic investigation. To bridge this gap, we introduce a data-driven framework that connects fish body shape dynamics to hydro-dynamic performance through compact morphospace parameterization and reduced-order modeling. Using CFD simulations of 15 fish species from the Digital Life Project database (www.digitallife3d.org/3d-model), we generate hydrodynamic datasets capturing the shape-flow relationship. Principal Component Analysis (PCA) extracts four dominant shape parameters from dorsal body profiles, which are then integrated into an Inverse-Design with Dynamic Mode Decomposition (ID-DMD) framework to model the resulting fluid dynamics. The resulting modal analysis suggests that locomotion strategies emerge from specific shape-flow interactions. We further demonstrate the frameworks utility through single- and multi-objective shape optimization, showing how it enables efficient exploration of the morphology-hydrodynamics relationship. This approach offers a novel analysis and design tool for understanding how biological form and motion interact with fluid mechanics, with applications ranging from bio-inspired vehicle development to evolutionary biomechanics.

Aortic stiffness is both heterogenous and anisotropic. Current non-invasive methods to estimate aortic stiffness are limited to characterizing the aortic tissue as isotropic due to the lack the techniques required to extract multi-axial strain from 3D dynamic images. Vascular deformation mapping (VDM) is a nonrigid image registration technique which has thus far been applied to map aortic growth using longitudinal imaging. In this study, we propose to use VDM to assess 3D aortic deformation by mapping diastolic and systolic images. During image registration process, penalty parameters are employed to fine-tune image alignment and penalize non-physiological deformations. These penalty parameters must be calibrated to ensure that VDM successfully reproduces multi-axial aortic motion patterns in health and disease. In this paper, we developed a calibration pipeline for these parameters using synthetic data. A rotation-free shell model was used to generate physics-based synthetic data on aortic motion incorporating patient-specific geometries, root motion, and blood pressure from a cohort of 14 subjects (healthy, Marfans syndrome and thoracic aortic aneurysm). An error metric was defined to quantify the quality of the VDM results. Furthermore, a k-means clustering technique was used to categorize the subjects into three clusters based on ascending aortic motion. Optimal penalty parameters were identified for each of the three clusters. The results indicated that patient clusters with smaller aortic root motion required larger rigidity penalty values. The calibrated parameters successively reduced errors in 3D displacement and multi-axial stretch compared to un-optimized VDM predictions, enhancing the accuracy of capturing aortic deformation from dynamic images. Among the different aortic regions, the ascending thoracic aorta exhibits the largest error reduction.

Patellar osteochondral allograft (OCA) transplantation is widely used to treat large full-thickness cartilage defects, yet long-term failure and reoperation rates remain high. Although surface congruity and osseous integration are emphasized clinically, cartilage thickness and mechanical compatibility between donor and recipient are not considered. Our previous work suggests that cartilage thickness mismatch can amplify local deformation at the graft boundary, potentially compromising graft longevity. This study investigates how combined mismatches in cartilage thickness and mechanical properties influence the local strain environment at the patellar OCA interface. Simplified two-dimensional axisymmetric finite element models of patellar OCA repair were developed in ABAQUS. Donor-to-recipient cartilage thickness ratios ranging from 0.33 to 3.25 were evaluated together with donor-recipient Youngs modulus mismatches (2.5-7.0 MPa). Cartilage was modeled using homogeneous linear elastic and functionally graded material formulations to account for depth-dependent stiffness. A compressive pressure of 1.0 MPa was applied to represent patellofemoral joint loading, and peak compressive and shear strains were quantified at the graft boundary. Cartilage thickness mismatch produced localized high-strain regions (HSR) of compressive and shear strain at the donor-recipient interface that were absent in thickness-matched constructs. Strain amplification increased with both thickness and mechanical property mismatch. Compressive strain exhibited directional asymmetry, with donor-side-thicker configurations producing greater amplification than recipient-side-thicker configurations. Incorporating depth-dependent cartilage stiffness reduced peak strain magnitudes but did not eliminate mismatch-driven strain amplification. These findings demonstrate that cartilage thickness and mechanical disparity can create HSR at the patellar OCA graft boundary that may predispose grafts to impaired integration and long-term failure.

Purpose: To identify, through iterative user-centered design, the auditory biofeedback requirements and sound preferences supporting gait training in children with cerebral palsy (CP), and to determine which feedback variables, sound mappings, and sound types yield clinically viable and movement-interpretable paradigms. Methods: The iterative process spanned two prototype phases. Prototype A comprised seven paradigms demonstrated to two experienced physiotherapists (Workshop 1A). Two of these were subsequently discarded owing to poor sound-movement interpretability and two were modified. Six paradigms were added to Prototype B, demonstrated to four children, five parents, and one therapist (Workshop 1B) and two therapists (Workshop 2B). Data were analyzed using systematic text condensation. Results: Within-child sound preferences varied with energy level and sensory state on a given day. Sound-movement interpretability tended to suffer for paradigms with greater acoustic complexity (e.g. computer-generated music). Therapists endorsed a repertoire spanning both movement quality and movement quantity targets. Participants independently proposed paradigms rewarding restrained and controlled movement, a feedback category absent from the current prototype. Conclusions: Session-level calibration is preferable to fixed sound profiles, requiring real-time interface support for paradigm adjustment. Acoustic complexity must remain subordinate to movement-sound interpretability. Paradigms targeting movement restraint are a development priority unaddressed in the literature.

Recent advancements in wearable sensors and machine learning show promise for estimating lower-body joint torques outside of laboratory settings. Inertial Measurement Units combined with Convolutional Neural Networks have proven effective for this task. However, the impact of different input data types and formats remains underexplored. This study investigates how variations in input data influence the prediction of lower-body joint torques during walking. Results indicate that while dataset choice causes only minor differences in prediction performance, the overall quality of the dataset plays a more critical role than the specific input variables in achieving accurate torque predictions using wearable sensors.

BackgroundPain with movement is common in adults with knee osteoarthritis (OA), but the effect of movement-evoked pain on gait is not well understood. This relationship is vital to understand as gait mechanics are associated with OA initiation and progression. Our current understanding of acute changes in pain and gait stems from extended bouts of walking, however these bouts likely dont represent real-world behavior. Therefore, understanding how gait changes with shorter, more intense bouts of activity may provide valuable insight into the pain experience. MethodsAdults with (n=19) and without (n=19) knee OA wore inertial measurement units (IMUs) while completing bouts of walking before and after two bouts of stair navigation (two flights). We tested whether pain and gait (speed, stride length, and lower extremity joint ranges of motion (ROM)) changed differently between adults with and without knee OA in response to multiple bouts of stair activity. FindingsThere were no significant interactions between group and stair bouts for any variable. When stratifying the OA group by those who did and did not experience pain, those who experienced a change in pain also had a greater change in early stance knee ROM in response to bouts of stairs. InterpretationThe observed changes suggest that knee kinematics may be more sensitive to acute changes in pain than gait speed or stride length. These differences were detectable using IMUs and therefore our results support the use of IMUs to measure concurrent pain and gait mechanics in less controlled and real-world settings.

Abstract Purpose: MyotonPRO (MTP) and time-harmonic elastography (THE) are increasingly used to assess muscle mechanical properties, yet they operate on fundamentally different physical principles. MTP measures composite MTP stiffness (N/m) through surface oscillations, while THE quantifies intrinsic shear modulus (THE stiffness, kPa) via propagating shear waves. This study aimed at systematically compare MTP and THE measurements in the vastus lateralis muscle across different contraction intensities and examine how the skin layer and subcutaneous fat (SLSF) thickness influence their relationship. Methods: Twenty-six healthy adults (15 males, 11 females; age 25 [SD 4] years) underwent MTP and THE measurements of the vastus lateralis at rest and during isometric contractions at 15% and 30% maximal voluntary contraction (MVC). Effects of contraction intensities on tissue properties were assessed using univariate analyses of variance with repeated measures. Associations between the different outcomes of THE and MTP technologies were explored using Pearson's correlations and partial correlation coefficients separately for each contraction intensity with adjustment of the SLSF thickness of participants. Results: Both technologies detected contraction intensity-dependent stiffening across all outcomes (p < 0.001). THE stiffness increased from 5.3 [1.2] kPa at rest to 15.6 [6.1] kPa at 30% MVC; THE wave attenuation increased from 0.83 [0.19] to 1.42 [0.36] s/m while MTP stiffness increased from 337.3 [49.3] N/m at rest to 529.4 [160.7] N/m at 30% MVC. Correlations between modalities were weak and condition-dependent. THE wave attenuation did not significantly correlate with any MTP outcome across conditions. Conclusion: MTP and THE detect contraction-induced stiffening through fundamentally different physical mechanisms and should not be regarded as interchangeable. Their correlation is modest at rest and breaks down (or reverses) during active contraction, with subcutaneous fat as a key modifying factor. Clinical trial number: Not applicable.

Hemodynamic forces play a key role in early cardiac morphogenesis, yet many computational studies assume Newtonian blood behavior. Here, we evaluate the impact of nonNewtonian shearthinning rheology on flow patterns, pressure distributions, and wall shear stress (WSS) during cardiac looping using idealized threedimensional models of the embryonic heart tube. Five geometries representing progressive looping stages, from a linear tube to an Sshaped configuration with ventricular ballooning, were analyzed under pulsatile flow using both Newtonian and powerlaw viscosity models. Across all stages, Reynolds numbers (Re {approx} 1-7) and Womersley numbers (Wo {approx} 0.3) indicated laminar, quasisteady flow consistent with embryonic conditions. Incorporating shearthinning rheology produced substantial deviations from Newtonian predictions, with peak systolic WSS differing by up to [~]40% and pressure drops by up to [~]20%. These effects were most pronounced in regions of increased curvature and geometric complexity. These findings demonstrate that nonNewtonian rheology significantly influences predicted hemodynamic environments during cardiac looping and should be incorporated into computational models aimed at understanding mechanobiological regulation of early heart development.